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Abstract 

In the quantization of simple cosmological models (minisuperspace models) described by the 
Wheeler-DeWitt equation, an important step is the construction, from the wave function, of a 
probability distribution answering various questions of physical interest, such as the probability 
of the system entering a given region of configuration space at any stage in its entire history. A 
standard but heuristic procedure is to use the flux of (components of) the wave function in a 
WKB approximation. This gives sensible semiclassical results but lacks an underlying operator 
formalism. In this paper, we address the issue of constructing probability distributions linked to 
the Wheeler-DeWitt equation using the decoherent histories approach to quantum theory. The 
key step is the construction of class operators characterizing questions of physical interest. Taking 
advantage of a recent decoherent histories analysis of the arrival time problem in non-relativistic 
quantum mechanics, we show that the appropriate class operators in quantum cosmology are 
readily constructed using a complex potential. The class operator for not entering a region of 
configuration space is given by the S-matrix for scattering oS a complex potential localized in that 
region. We thus derive the class operators for entering one or more regions in configuration space. 
The class operators commute with the Hamiltonian, have a sensible classical limit and are closely 
related to an intersection number operator. The definitions of class operators given here handles 
the key case in which the underlying classical system has multiple crossings of the boundaries of 
the regions of interest. We show that oscillatory WKB solutions to the Wheeler-DeWitt equation 
give approximate decoherence of histories, as do superpositions of WKB solutions, as long as the 
regions of configuration space are sufficiently large. The corresponding probabilities coincide, in 
a semiclassical approximation, with standard heuristic procedures. In brief, we exhibit the well- 
defined operator formalism underlying the usual heuristic interpretational methods in quantum 
cosmology. 

PACS numbers: 04.60.-m, 04.60.Kz, 04.60.Ds, 03.65.Yz 
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I. INTRODUCTION 



A. Opening Remarks 

A question in quantum gravity that continues to attract considerable interest is the 
problem of quantizing and interpreting the Wheeler- De Witt equation of minisuperspace 
quantum cosmology [1-5], 



Although the setting of this problem is simple cosmological models with just a handful 
of homogeneous parameters, the techniques employed in answering this question may be 
relevant to general approaches to quantum gravity, such as the loop variables approach or 
causal set approach. This is because the central difficulty in consistently quantizing and 
interpreting the Wheeler-DeWitt equation is the absence of a variable to play the role of 
time and all approaches to quantum gravity must confront this issue at some stage [6]. 

A frequently studied example consists of a closed FRW cosmology with scale factor a = 
and a homogeneous scalar field (f) with (infiationary) potential [7]. The Wheeler-DeWitt 
equation for this model is 



Given suitable boundary conditions, one can solve this equation for the wave function "if {a, 0) 
and attempt to use it to answer a number of interesting cosmological questions. There 
are many such questions: Is there a regime in which the universe behaves approximately 
classically? What is the probability that the universe expands beyond a given size Oq? What 
is the probability that the universe has a certain energy density at a given value of the scale 
factor? What is the probability that the universe's history passes through a given region A 
of configuration space, characterized by certain ranges of a and 0? 

Such questions, of necessity, do not involve the specification of an external time. Classi- 
cally, the system's trajectories in minisuperspace may be written as paths but 
here t is a convenient but unphysical parameter that labels the points along the paths - it 
does not correspond to the physical time measured by an external clock. The absence of a 
physical time is refiected in the quantum theory by the fact that the quantum state obeys 
the Wheeler-DeWitt equation Eq.(l.l), not a Schrodinger equation, and it is this difference 
that presents such a challenge to conventional quantization methods. 





(1.2) 
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Although very plausible heuristic semiclassical methods exist for formulating and answer- 
ing the above questions (in particular, the WKB interpretation [1]), it is of interest to see 
whether or not there exists a precise and well-defined quantum-mechanical scheme underly- 
ing these heuristic methods. We are not looking for high standards of mathematical rigour 
- just a standard quantum-mechanical framework of operators, inner product structures 
etc, obeying some reasonable requirements. The purpose of this paper is to show, building 
on earlier attempts [8-15], that such a framework exists in the context of the decoherent 
histories approach to quantum theory. 

B. Inner Products and Operators for the Wheeler-DeWitt Equation 

In the Wheeler-DeWitt equation (1.1), H is the Hamiltonian operator of a minisuperspace 
model with n coordinates g", and is typically of the form 

H = -V^ + U{q) (1.3) 

where V is the Laplacian in the minisuperspace metric /„;,, which has signature (— • • ). 
It is naturally linked to the current, 

Ja = i (^*^a* - (1.4) 

which is conserved 

V-J = (1.5) 
Closely associated is the Klein-Gordon inner product defined on a surface E 

(*, ^)kg = '1^^" {^*~da^ - ^* (1.6) 

where da"' is a surface element normal to S. In fiat space with a constant potential, the 
Whcclcr-DeWitt equation is just the Klein-Gordon equation. Its solutions may be sorted 
out into positive and negative frequency in the usual way. With a little attention to sign, it 
is then possible to use components of the current to define probabilities. 

However, in general, it is not possible to sort the solutions to the Wheeler-DeWitt equa- 
tion into positive and negative frequency. This is one manifestation of the problem of time 
and more elaborate methods are required to associate probabilities with the Wheeler-DeWitt 
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equation. There are two main issues: finding an inner product, and then finding suitable 
operators. 

The issue of finding a suitable positive inner product is reasonably straightforward and 
goes by the name of Rieffel induction, or the induced (or physical) inner product [9, 16]. We 
consider first the usual Schrodinger inner product, 

(^i|^2) = J dy^l{q)^2{q) (1.7) 
We then consider eigenvalues of the Wheeler-DeWitt operator 

i/I^Afe) = AI^Afe) (1.8) 

where k labels the degeneracy. These eigenstates will satisfy 

{^yk'\'^Xk)^S{X-X')S{k-k') (1.9) 

from which it is clear that this inner product diverges when A = A'. The induced inner 
product on a set of eigenstates of fixed A is defined, loosely speaking, by discarding the 
5-function 5{X — A'). That is, the induced (or physical) inner product is then defined by 

This procedure can be defined quite rigorously, and has been discussed at some length in 
Refs.[9, 16]. It is readily shown that the induced inner product coincides with the Klein- 
Gordon inner product when a division into positive and negative frequencies is possible, 
with the signs adjusted to make it positive. (This is described in Appendix A). 

Turning now to the construction of interesting operators, the interesting dynamical vari- 
ables associated with the Wheeler-DeWitt equation are those that commute with H. This 
is because the constraint equation is related to reparametrization invariance - which is re- 
flected in the absence of a physical time variable - and we seek operators which are invariant. 
A wide class of operators commuting with H are of the form 



A= dt Bit) (1.11) 

J —oo 

which clearly commute with H since 

/oo 
dt B{t + s) 
-oo 

= A (1.12) 



Many examples are given in Refs.[5, 11, 16]. However, another way of constructing such 
operators involves taking products, 

oo 

H B{t) (1.13) 

t=— oo 

which may be shown to commute with H using essentially the same argument [12], but clearly 
further mathematical detail is required to give meaning to the infinite product. (Here t is 
the unphysical parameter time). 

Given these prescriptions for inner products and operators, one may then attempt to 
construct operators and probabilities implementing some of the questions mentioned above. 
We will focus on the following general question: Given a solution ^ to the Wheeler-DeWitt 
equation, what is the probability of finding the system in a region A of configuration space, 
or of crossing a surface E, at any stage in the system's history? The question is depicted in 
Figure 1. 

Questions involving surface crossings are not unlike more familiar questions in non- 
relativistic quantum mechanics, where there is a physical time parameter, but the key 
difference in quantum cosmology is that even classically, a given trajectory will typically 
cross a fixed surface more than once. It is precisely these sorts of issues that need to be 
carefully phrased in the quantum theory. Whilst the operator formalism briefly outlined 
above has been used to address such questions [5], the problems of characterizing properties 
of trajectories and surface crossings in quantum theory is naturally accommodated in the 
decoherent histories approach to quantum theory. 

C. The Decoherent Histories Approach 

A more general approach to quantizing and interpreting the Wheeler-DeWitt equation 
is provided by the decoherent histories approach to quantum theory, and this approach will 
be the focus of this paper [17-25]. In this approach, probabilities are assigned to histories 
using the formula 

p{a)^TT{C^pCi) (1.14) 

Here, p is the initial state (in our case a pure state) and Ca is a class operator characterizing 
the histories a of interest. In non-relativistic quantum mechanics these class operators are 
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given by time-ordered strings of projection operators, 

{t2)Pa,{tl) (1-15) 

(or by sums of such strings) . The class operators alway satisfies the condition 

J2Ca=l (1.16) 

a 

For the reparamterization invariant systems considered here, the definition of the class op- 
erators is more subtle, and wc return to this below. 

Because of interference between pairs of histories, probabilities cannot always be assigned. 
To check for interference we therefore consider the decoherence functional, 

D{a,a')^Tr(c^pCl') (1.17) 

When 

D{a,a')^0 (1.18) 

for all pairs of histories in the set with a ^ a', we say that there is decoherence of the set 
of histories and probabilities may be assigned Eq.(1.14). When there is decoherence, it is 
easily seen from Eq.(1.16) that the probabilities Eq.(1.14) are also given by 

p{a) = Tr{C^p) = Tr{Cip) (1.19) 

Decoherence guarantees that this expression is real and positive, even though the class 
operators are not positive or hermitian operators in general. 

The structure of the dccoherent histories approach is very general and may be applied 
to a wide variety of situations, given an initial state, class operators, and a suitable inner 
product structure. For the apphcation to the Wheeler- De Witt system considered here, the 
initial state is take to be a solution to the Wheeler-DeWitt equation and the inner product 
is the induced inner product described above. 

The most crucial element is the specification of the class operators Ca- In non-relativistic 
quatum theory, the class operator is generally a string of projection operators like Eq.(1.15), 
or a sum of strings of such operators, but this basic structure can be generalized in various 
ways. The product in the string can be taken to be continuous time [25]. Also, it is 
often valuable, sometimes essential, to allow the projectors to be replaced by more general 
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operators, such as POVMs. The class operators must also properly characterize the histories 
one is interested in. It is not always obvious how to do this but useful clues often come from 
looking at the classical analogue of the class operator (where all the projectors commute). 

Here, we are interested in histories which enter a region A of configuration space, or 
which cross a surface S, but without regard to time. This absence of a physical time 
variable seems particularly challenging, given that time seems to be central to the definition 
of non-relativistic analogue, Eq.(1.15). Closely related to this is the role of the constraint 
equation, Eq.(l.l). As noted already, these two features are directly related to the underlying 
symmetry of the theory - reparametrization invariance - and this symmetry is respected if 
the class operators commute with the constraint, 

[H,C^]^0 (1.20) 

Eq.(1.20) is in keeping with standard procedures of Dirac quantization [5]) and also ensures 
that the class operators have a sensible classical limit [12]. 

Some partially successful attempts to construct such class operators have been given pre- 
viously [9-12], but ran into various problems to do with the Zeno effect and with compati- 
bility with the constraint equation. The main aim of this paper is to give fully satisfactory 
definitions of class operators for quantum cosmological models and explore their decoherence 
properties and probabilities. 

D. Some Properties of the Wheeler-DeWitt Equation and the WKB Interpreta- 
tion 

To prepare the way for the full decoherent histories analysis of quantum cosmology, it 
is important to discuss some properties of the Wheeler-DeWitt equation and review the 
commonly used heuristic semiclassical interpretation of the wave function, since a proper 
quantization must recover this structure in some limit. 

The Wheeler-DeWitt equation Eq.(1.3) is a Klein-Gordon equation in a curved configura- 
tion space with indefinite metric fab{(l) and potential U{q) which can be positive or negative. 
The curvature effects of the metric are not significant in relation to the issues addressed in 
this paper, so we will assume for simplicity that the metric is fiat. 

The classical constraint equation corresponding to the Wheeler-DeWitt equation Eq.(1.3) 
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IS 



(1.21) 



from which one can see that the classical trajectories are timelike in the region U > 
and spacelike in U < 0. (The timelike direction is that of increasing scale factor and the 
spacelike directions correspond to matter and anisotropic modes). The quantum case has 
analogous features. In simple models such as Eq.(1.2), the character of the solutions to the 
Wheeler- De Witt equation depends on the sign of U. For large scale factors, U > and 
the wave function is oscillatory, corresponding, very loosely, to a quasiclassical regime, and 
for small scale factors, U < 0, and the wave function is exponential, corresponding to a 
classical forbidden regime. However, there are certain types of models (such as those with 
an exponential potential for the scalar field), in which the identification of the oscillatory 
and exponential regions depends also on whether the constant U surfaces are spacelike or 
timelike [26] . We will not address this here. 

One can also associate a Feynman propagator Gp with the Wheeler-DeWitt operator. 



where e — > 0+. (Numerous propagator- like objects of this type have been considered in 
quantum cosmology [27]). Locally, on scales smaller than the scale on which the potential 
U significantly varies, the propagator Gpiq, q') will be essentially identical in its properties 
with the analogous object for the relativistic particle for fiat space. Therefore, in the region 
U > 0, for points q and q' which are timelike separated, it will propagate positive frequency 
solutions to the future and negative frequency solutions to the past. It will be exponentially 
suppressed for initial and final points q and q' which are spacelike separated (this is the 
familiar "propagation outside the lightcone" effect). For [/ < 0, similar statements hold but 
with timelike and spacelike reversed. Similar features hold on larger scales in a semiclassical 
approximation. These properties are important to understand the class operator constructed 
below. 

Very plausible but heuristic answers to questions concerning crossing surfaces and entering 
regions are readily found using the WKB approximate solutions to the Wheeler-DeWitt 
equation and the Klein-Gordon current [1]. In the oscillatory regime, the Wheeler-DeWitt 




'OO 



dt e 



—I 



(1.22) 



H-it 
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equation may be solved using the WKB ansatz 

* = Re'^ (1.23) 

where the rapidly varying phase S obeys the Hamilton-Jacobi equation 

{VSf + [/ = (1.24) 

and the slowly varying prefactor it! obeys 

V- (li^l^V^) = (1.25) 

The latter equation is just current conservation for the WKB current, 

J=\R\^VS (1.26) 

Wave functions of the WKB form Eq.(1.23) indicate a correlation between position and 
momenta of the form 

p = VS (1.27) 

and this suggests that the wave function Eq.(1.23) corresponds to an ensemble of classical 
trajectories satisfying Eq.(1.27). The current J may then be used to define a measure on 
this set of trajectories. For example, we consider a surface E and choose the normal to 
the surface so that n ■ VS > 0. Then the probability of the system crossing the surface is 

taken to be 

p(E) = j r-'q n'^Ja (1.28) 

Of particular interest here is the probability of entering a region A. This is clearly related 
to the flux at the boundary of the region. The current will typically intersect the boundary 

of A twice. However, wc can split the boundary E of A into two sections: Ei„ at which the 
current is ingoing and Eq^^ at which the current is outgoing. The probability of entering A 
is may then be expressed in the two forms 

P(A) = - I r-'q n^Ja 

(T-^q rfJa (1.29) 



out 



where here we have defined the normal to point outwards. See Figure 2. These two forms 
are equivalent since the current is locally conserved. 
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Typically, little is said of the regions in which the wave function is exponential, except 
that they are similar to tunneling regions in non-relativistic quantum mechanics, so are 
classically forbidden in some sense. However, the wave function is not necessarily small in 
these regions, so there is surely more to it than this. In Ref.[28], it was noted that, unlike the 
oscillatory regions, the exponential regions do not indicate a correlation between position 
and momenta of the form Eq.(1.27), and it seems that this should be significant in some 
way. 

Whilst the WKB interpretation is very plausible and adequate for most situations of 
interest, it leaves many questions unanswered. The main issue is to understand the op- 
erator origin of these probabilities, in terms of the language of operators commuting with 
H developed above. Furthermore, what can one say about superpositions of WKB states? 
Can interference between them be neglected? Also, what can one say about the expo- 
nential, classically forbidden regions? Is there a more precise way of saying that they are 
non-classical? 

E. This Paper 

The purpose of this paper is to present a decoherent histories quantization of the Wheeler- 
DeWitt equation, and in particular to exhibit exactly defined class operators which charac- 
terize histories entering regions of configuration space or crossing surfaces, without reference 
to an external time. The key idea is to use a complex potential to define the class operator 
for not entering a region A in configuration space. In particular, we take the class operator 
for not entering to be the S matrix describing scattering off a complex potential localized 
in A. This turns out to have all the right properties and to give physically sensible results 
for decoherence and probabilities. 

In Section 2, to explain and motivate the use of a complex potential we review the 
use of such potentials in the decoherent histories analysis of the arrival time problem in 
non-relativistic quantum mechanics. In Section 3, we describe the construction of class 
operators for the Wheeler- De Witt equation using a complex potential. The properties of 
the class operator for entering a region A arc described in Section 4. The class operator 
has a sensible classical limit which, crucially, registers just one intersection of a trajectory 
with a surface, even when the trajectory intersects twice or more. In the quantum case, the 
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class operator is an operator describing the ingoing flux across the boundary of the region, 
an expected result on semiclassical grounds, and is closely related to an intersection number 
operator. 

In Section 5, some simple one and two dimensional examples are considered in detail to 
confirm the properties of the class operator outlined in Section 4. We also confirm that the 
formalism gives sensible and expected results for surface crossings of the relativistic particle 
in flat space. In Section 6 we consider the WKB regime and show that the decoherent 
histories analysis reproduces the expected heuristic interpretation of the wave function. We 
summarize and conclude in Section 7. Some properties of the inner product structure of the 
Klein-Gordon equation are summarized in Appendix A. 

II. THE ARRIVAL TIME PROBLEM IN NON-RELATIVISTIC QUANTUM THE- 
ORY 

In this section we summarize some of the key features of the arrival time problem in 
non- relativistic quantum mechanics [29]. These details are very relevant to the quantum 
cosmology case and in particular motivate the use of complex potentials in the definition of 

class operators. 

In the onc-dimensional statement of the arrival time problem, one considers an initial 
wave function 1-0) concentrated in the region x > and consisting entirely of negative 
momenta. The question is then to find the probability Il{T)dT that the particle crosses 
X — between time r and T + dr. See Figure 3. The canonical answer is the current density 



where Itpr) is the freely evolved state. This has the correct semiclassical hmit, but can be 
negative for certain types of states consisting of superpositions of different momenta (back- 
flow states). It is of interest to explore whether this simple semiclassical result emerges from 
more elaborate measurement-based or axiomatic schemes. Many such schemes naturally 
lead to the use of a complex potential —iV{x) in the Schrodinger equation [30], where 




(2.1) 



V{x) = Vo9{-x) 



(2.2) 



With such a potential, the state at time r is 



V^(t)) = exp {-iHr - Voe{-x)T) 



(2.3) 
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where H is the free Hamiltonian. The idea here is that the part of the wave packet that 
reaches the origin during the time interval [0, r] should be absorbed, so that (■0(t)|'0(t)) 
is the probability of not crossing x — Q during the time interval [0,r]. The probability of 

crossing between r and t + dr is then 

n(r) = -^(^(r)|^(r)) (2.4) 

This may be approximately evaluated in the limit Vq <^ E (where E is the typical energy 
scale) , which is the limit of neghgible reflection off the complex potential. The result is 

n(T) = 2K, f dt e-2^°(^-*) J{t) (2.5) 
Jo 

The probability for crossing during a finite interval [ti, T2] is then given by 

P(r2,ri)= f'drUir) (2.6) 



and if this time interval is sufficiently large compared to I/Vq, we have 

p(t2,ti) ^ f dt J{t) (2.7) 

which means that the dependence on the potential drops out entirely at sufficiently coarse 
grained scales, and there is agreement with semiclassical expectations involving the current 
J{t). 

In the decoherent histories analysis of the arrival time problem [31, 32], we first consider 
the construction of the class operator C„c for not crossing the origin during the finite time 
interval [0, r]. We split the time interval into N parts of size e, and the class operator is 
provisionally defined by 

Cnc = Pe-'"'P- P (2.8) 

where there are A'" + l projections P — 9{x) onto the positive x-axis and N unitary evolution 
operators in between. One might be tempted to take the limit — > 00 and e — > 0, but this 
yields physically unreasonable results. This hmit actually yields the restricted propagator 
in a: > 0, 

Cnc = 9'r(T,0) 

= Pe-^P"P^ (2.9) 
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This object is also given by the path integral expression 

(a;i|^,(r,0)|a;o) = / Vxexp{iS) (2.10) 

J r 

where the integral is over all paths from a;(0) = xq to x{t) = xi that always remain in 
x{t) > 0. However, the class operator defined by Eq.(2.9) has a problem with the Zeno 
effect - it consists of continual projections onto the region a; > and as a result the wave 
function never leaves the region. This is reflected in the fact that the restricted propagator 
Qr is unitary in the Hilbert space of states with support only in x > 0. This is a serious 
difficulty which has plagued a number of earlier works in this area [33-35] . 

To avoid the Zeno effect, the key is to keep e finite, and in particular, standard wisdom 
suggests that e > l/AH, where H is the free particle Hamiltonian. The class operator 
Eq.(2.8) is not easy to work with for finite e, but fortunately an extremely useful result 
of Echanobe et al. comes to the rescue [36]. They argued that the string of operators in 
Eq.(2.8) is in fact approximately equivalent to evolution in the presence of the complex 
potential —iV introduced above. That is. 



Pe-'"'P---e-''''P !^ ex.p{-iHoT-Voe{-x)T) (2.11) 

They argued that this is valid for AH <^ Vq and T^oe ^ 1, but there is evidence that this 
connection is valid more generally [37]. In any event it strongly suggests that the class 
operator C„c, normally defined by a string of projection operators, is justifiably defined 
instead using a complex potential. That is, we define 

Cnc = exp{-iHoT-Vo9{-x)T) (2.12) 

The subsequent decoherent histories analysis was described in detail in Ref.[31]. The corre- 
sponding class operator for crossing during a time interval [ti, T2] was shown to be 

Cc(t2, Ti) = r dt e-^^o(^-*)v^e-^^°*-^* (2.13) 

and again in the approximation Vq -C and for time intervals greater than 1/Vq this may 
be shown to have the simple and appealing form 

Cc(r2,ri) ^ J^^'dt e^'""'^ (P 5ixt) + 5{xt) p) (2.14) 

which is now independent of the Vq. (These definitions of class operators differ by a factor 
of exp(— iifr) from those defined in Section IC, but this does not make any difference in the 
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decoherence functional and probabilities). With this class operator, one can show that there 
is decoherence of histories for a variety of interesting initial states, such as wave packets (but 
not for superposition states with backflow), and for such states, the general result Eq.(1.19) 
implies that the probability for crossing is simply 



which agrees precisely with the semiclassically expected result, Eq.(2.7). 

In summary, the decoherent histories analysis of the arrival time problem in non- 

relativistic quantum mechanics indicates that it is reasonable to define class operators for 
not entering a space time region using a complex potential, as in Eq.(2.12), and that such 
a definition gives sensible a semiclassical limit, independent of the potential, at sufficiently 
coarse grained scales. 

III. CONSTRUCTION OF THE CLASS OPERATORS USING A COMPLEX PO- 
TENTIAL 

We now come to the central issue concerning this paper, which is the construction of class 
operators for the decoherent histories analysis of the Wheeler-DeWitt equation. 

A. Class Operators for a Single Region 

We seek class operators for the system Eq.(1.3) describing histories which enter or do not 
enter the region A, without specification of the time at which they enter. It is easiest to 
first focus on the class operator Ca for not entering and the class operator for entering is 
then given by 



The earhest attempts to define class operators for the Wheeler-DeWitt equation involved 
defining C/^ as a sort of propagator obtained by integrating restricted propagators of the 
form Eq.(2.10) over an infinite range t (now regarded as the unphysical parameter time) 
[9-11]. However, in addition to having problems with the Zeno effect (which was not in 
fact appreciated in these earlier works) , such constructions are difficult to reconcile with the 




(2.15) 



Ca — 1 — Ca 



(3.1) 
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constraint equation and some ad hoc modifications of the basic construction were required 
to give sensible answers. 

A rather different approach to constructing (7a was given in Ref.[12]. This was again 
problematic, but we review the construction here, since it is readily modified to yield a 
successful definition of the class operators. We denote by P the projector onto A and P 
the projector onto the outside of A. Our provisional proposal for the class operator for 
trajectories not entering A is the time-ordered infinite product, 

oo 

= n ^(^) (3-2) 

t=—oo 

where t is the unphysical parameter time. Subject to a more precise definition, given shortly, 
this object has the required properties. It is a string of projectors. Classically, it is equal 
to one for trajectories which remain outside A at every moment of parameter time. Also, it 
commutes with if, at least formally. 

To define this more precisely, we first consider the product of projectors at a discrete set 
of times, ti, ti + e, ti + 2e, ■ ■ ■ ,ti + ne = t2- We define the intermediate quantity, CA{t2, ti) 
as the continuum limit of the product of projectors, 

CAih, h) = lim P{h)P{h - e) ■ ■ ■ P{ti + e)P(ti) (3.3) 

where the limit is n — > oo, e — > with t2 — ti fixed. The desired class operator is then 

Ca= lim CA{t2,ti) (3.4) 

t2 — ^0O,tl — > — 00 

The class operator operator is clearly closely related to the restricted propagator g^. (the 
generahzation of Eq.(2.10)) in the region outside A, since we have 

C4t2, h) = e'"'' gr{t2, h) e-^^*i (3.5) 

and therefore 

(5a = lim e'^^^ gr{t2,ti) e-'"^' (3.6) 
The class operator commutes with H. This is because, from Eq.(3.5) 

e''''CAit2, t{)e-'"' = e'^(*^+^) ^,(^2, h) e-'"^''+'^ (3.7) 

This becomes independent of s SkS t2 ^ oo, ti — > — oo, hence 

[i/,C'A]=0 (3.8) 
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The problem with this definition, however, is that it suffers from the Zeno effect, exactly 
like the analogous expression Eq.(2.8) (in the limit e — > 0) in the non-relativistic arrival time 
problem. But the key idea here is that we may also get around the problem in the same 
way, using a complex potential. That is, we "soften" the restricted propagator and make 
the replacement 

gr{t2: ti) ^ exp {-iH{t2 - U) - V{t2 - h)) (3.9) 

where V{q) = Vof^{q). Here, Vq > is a constant and /a(9) is the characteristic function 
of A, so is 1 in A and outside it. Or we may equivalently write V = VqP, where recall, P 
is the projector onto A. This means that our new definition for the class operator is 

Ca= lim e*^*^ exp {-i{H - iV){t2 - h)) e"'^*! (3.10) 

The class operator thus defined has an appealing form: it is the S-matrix for system with 
Hamiltonian Eq.(1.3) scattering off the complex potential V. As required, it commutes with 
H. This is the most important definition of the paper. 

Now an important observation. The class operator derived above has been defined as 
a time-ordered product of operators in which the direction of parameter time increases 
from right to left. However, since parameter time is unphysical, there is absolutely no 
reason why the parametrization should not run in the opposite direction. This produces an 
operator which is the hcrmitian conjugate of Eq.(3.10). As one can see from Eq.(1.19), this 
makes no difference in the final expressions for probabilities. In fact, it is most natural to 
define the class operator in such a way that it is invariant under reversing the direction of 
parametrization. We thus define a modified class operator which is hermitian: 

C'A^l{CA + Ci) (3.11) 

In what follows, for simplicity, we will primarily work with the non-hermitian class operator 
Eq.(3.10) and revert to the hermitian one, Eq.(3.11) where appropriate. The difference 
between them will turn out to be significant only for the class operator for two or more 
regions. 

We now cast the above class operator in a more useable form. The following identities 
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are readily derived: 

r't2 



rt-z 

^-i{H-iV)(t2-ti) ^ ^-iH{t2-ti) _ / ^-iH(t2-t)Y^-i{H-iV)(t-ti) (3.12) 

iH{t2-h) _ ^-^(H-iV){t2-t)y^-^H{t-H) ^g^^^g) 

Ju 



— e 



Inserting the second expression in the first, we obtain 

g-j(if-iy)(t2-tl) ^ ^-iH{t2-tl) _ ^-iH{t2-t)Y^-iH{t-tl) 

Jti 

+ r dt rc^se-'^(*^-*Ve-^(^-^^)(*-^Ve-^^(^-*i) (3.14) 
Inserting in the expression for the class operator (3.10) and taking the limit, we obtain 

/oo 
dt V{t) 
-OO 

/OC /'t 
dt / ds e^^Ve-'(^-'^)(*-^Ve-'^^ (3.15) 
-oo J —oo 

The class operator for entering the region is therefore given by 

/cx) roc rt 

dt V{t) - dt ds e^^*T/e-^(^-^^)(*-^)ye-^^^ (3.16) 
-oo J — oo J — oo 

This is an exact and useful form for the class operator for entering A, and it is easily 
confirmed that it is of the form Eq.(l.ll) so commutes with H. 

Now we use a simple but useful semiclassical approximation. Noting that V — VqP, 
where P recall projects into A, note that the expression 

y^-iiH-iV){t-s)y (3_;^7) 

describes propagation with the complex Hamiltonian H — iV between two points that lie 
inside A. This can easily be represented by a path integral in which it seems plausible that 
the dominant paths between this two end-points will lie entirely inside A, as long as the 
boundary is reasonable smooth and the semiclassical paths are not too irregular. If this is 
true, we may replace V = VqP by the constant complex potential V = Vo, that is. 

Propagation with a complex potential in Eq.(3.17) will also involve reflection off the bound- 
ary of the region, in which case the semiclassical approximation Eq.(3.18) may fail. However, 
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reflection is smaU for sufficiently smaU Vq [30, 31], and we will see this in more detail in Sec- 
tion 5. Hence we expect te semiclassical approximation Eq.(3.18) to hold for small Vq. 

With this useful approximation (which does not affect the fact that the class operator 
commutes with H), we have 

/oo roo pt 

dt V{t) - dt ds y(i)y(s)e-^°(*-") (3.19) 
oo J —oo J —oo 

Again using V — VqP, this is easily rewritten 

/oo /»^ 
dt P{t) I ds Vbe-^°(*-'')p(s) (3.20) 
■oo J —oo 

This is the main result of the paper: a class operator commuting with i?, describing histories 
which enter the region A. It is valid in the approximation Eq.(3.18), which is sufficient to 
cover the key case of histories which, classically, intersect the boundary of A twice. 

For systems whose classical paths intersect the surfaces of interest more than two times, 
the semiclassical approximation Eq.(3.18) will not be valid, but the exact result Eq.(3.16) 
may still be used. It may also be of interest to explore higher order semiclassical approx- 
imations obtained by iterations of the basic result Eq.(3.14). This will not be explored 
here. 

B. Class Operators for Two Regions 

We now extend the above results to the case of histories which enter two regions, Ai, A2 
where we define Pi and P2 to be the projection operators onto these regions. We proceed 
as follows. The class operator Eq.(3.15) is defined for any complex potential V so we may 
use a potential 

y = 14 + T/2 

= FioPi + 1^20^2 (3.21) 

which has support only in the regions Ai, A2. Here Vio and V20 are positive constants. This 
potential inserted in Eq.(3.15) defines the class operator C'aiA2 for histories which remain 
outside of both Ai and A2. The negation of this class operator therefore describes histories 
which enter Ai or A2, or both, and we may write 

C'aiA2 = 1 ~ C'ai — — C'aiA2 (3.22) 
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Here, Cai and Caj are the class operators for entering the regions Ai, A2 and are given 
by expressions of the form Eq.(3.16). We may thus deduce the form of the class operator 
CaiA2 for entering both regions. 

The class operator Eq.(3.15) with the potential Eq.(3.21) is 

/oo 
dt {V,{t)+V2{t)) 
-CO 

/oo rt 
dt ds e'''\Vi + V2)e-'^''-'^^'+^'^^^'-'\Vi + V2)e-'''' (3.23) 
-00 J —00 

It seems reasonable to make the semiclassical approximations 

y^^.i(^H-i{V,+V^))(t-s)y^ _ y^^.i(^H-iV,){t-s)y^ ^334) 
y^^-i^H-i{V,+V^)){t-s)y^ _ y^^^i^H-iV^){t-s)y^ (g ^g) 

These approximations are very similar to Eq.(3.18). Eq.(3.24) is essentially the notion that 
the propagator between two points in Ai will involve negligible contribution from paths 
that enter A2. Similarly for Eq.(3.25). With these approximations we can easily identify 

the terms Cai and Caj and we deduce 

/oo pt 
dt / ds e'my^e-'^H-iiVMm-s)y^^-iHs 
-00 J —00 

/oo pt 
dt / ds e'Hty^e-'^H-i{V,+V,))(t-s)y^^-iHs (g ^g) 

-00 J —00 

To estimate the propagator 

y^f.-iiH-i{v,+V2)){t-s)y^ (3.27) 

we need a semiclassical approximation more elaborate than Eq.(3.18), since it involves prop- 
agation from an initial point in Ai with potential to a final point in A2 with potential V2 
but with zero complex potential outside these regions. The path decomposition expansion 
[38-40] is naturally adapted to this problem, since it is a tool for breaking up a propagator 
into its properties in different spatial regions. We will simply write down the result of using 
this expansion, with some brief justification. 

p^e-i{H-ny,+v,))(t-s)p^^ j' ds' j* dt' P2e-^^(*-*')P2e-^^(*'-^')Pie-'^(^'-^)Pi 

X exp(-V2o(t-0 -^io(s'-s)) (3.28) 
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This expression is actually exact when the complex potential is absent, as is easily verified. 
Since Pi and P2 represent the fluxes at the boundary of Ai and A2, one can easily identify the 
sections of the propagation which, semiclassically, are in Ai, outside Ai and A2, and in A2. 
The exponential suppression factors involving the complex potential for propagation inside 
Ai and A2 are inserted using a semiclassical approximation along the lines of Eq.(3.18). 
Inserting in Eq.(3.26), we thus obtain 

/oo pt pt pt 

dt ds ds' / dt' V2{t)P2{t')Pi{s')Vi{s) 

-00 J —00 Js J s' 

X eK^{-V2o{t-t')-ViQ{s' -s)) 

/OO nt rt pt 

dt ds ds' / dt' Vi{t)Pi{t')P2is')V2{s) 
-00 J —00 Js J s' 

X exp(-V^io(^-0-^2o(s'-s)) (3.29) 

IV. PROPERTIES OF THE CLASS OPERATORS 

We now examine the properties of the class operators Eq.(3.20), (3.29) and conflrm that 
they give the desired results. 

A. Single Region 

It is easy to show that the class operator for a single region Eq.(3.20) has a sensible 
classical hmit. Classically, P{t) is a function on classical trajectories with P{t) — 1 when 
the classical trajectory is in A and is zero otherwise. Suppose a given trajectory q^{t) enters 
A at some stage in its history so intersects the boundary twice (recalling we have essentially 
assumed no more than two intersections in the semiclassical approximation Eq.(3.18)). For 
a given choice of parametrization of the trajectory, it enters A at parameter time ta and 
leaves at time > ta (see Figure 4). (We may parameterize it in the opposite direction, 
with the same ultimate result). Then the derivative of P{s) is 

P{s)^5{s-ta)-5{s-t,) (4.1) 

and we have 

Ca= r dt I ds V^e-'"''^'-'^ [5{s - to) - 5{s - t^)] (4.2) 

Jta ''ta 
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Since s < t < U, the second 5-function in Eq.(4.2) makes no contribution to the integral. 
This is exactly the desired property - the expression for Ca registers only the first intersection 
of the trajectory with the boundary, but not the second intersection. The integral is easily 
evaluated with the result, 

Ca = 1 - 6-^°^'"-'^^ (4.3) 
This is approximately 1, as required, as long as 

Voih - ta) » 1 (4.4) 

We now give the broad picture in the quantum case and confirm some of the details in 
some simple models in the next section. Suppose we operate with Ca on an eigenstate I^'a) 
of H. The time integrals in Eq.(3.20) are easily carried out and the result is 

CaI^a) = 2TrVo5{H - A) PGyP|*A) (4.5) 



where 



G 



V 



POD 

/ dt e-*(«-^)*-^o* (4.6) 
Jo 

(4.7) 



H-X-iVn 



and we have used 



1 r°° 

S{H-X) = — / dt e-^(^-^)* (4.8) 
The P term is a current operator on the boundary E of A and may be written 

P = i[H,P] 

= -PnS^iq) - Ss{q)pn (4.9) 

where 

5^(q) = / (T-ig \q){q\ (4.10) 

is a 5-function operator in the surface S and Pn = n"'Pa is the component of the momentum 
operator normal to E with the outward pointing normal. The normal will depend 
on q in general so there may be an operator ordering issue in making pn hermitian. The 
difference between different ordering will involve a term of the form Van", essentially the 
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extrinsic curvature of E, and this will be small if A is sufficiently large and its boundary 
reasonably smooth. Note that it is sometimes also convenient to write P as 

P^i j^cT-^q \q)X{q\ (4.11) 

It is very useful to separate the current operator at the boundary into ingoing and out- 
going parts according to the sign of Pn at the boundary 

P = (P),„ - {PUt (4.12) 

where {P)in consists of the components of P with incoming momentum 

{P)in = -PnO{-Pn)SE{q) - S^{q)pn9{-pn) (4.13) 

and similarly for {P)out- Examples of these expressions in particular models will be given in 
the following sections. The restriction to positive or negative Pn means that it is generally 
difficult to express (P)i„ and {P)out in the form Eq.(4.11), involving the derivative pn = —idn, 
unless operating on a state (such as a WKB state) with simplifying properties. Note that 
these definitions require only that the local flux operator on a given surface can be split into 
ingoing and outgoing parts. To do so globally on a family of surfaces is generally impossible, 
essentially due to the problem of time, but fortunately this is not required here. 

The quantity Gy has the form of a Feynman propagator, Eq.(1.22), with Vq playing the 
role of the "ie prescription" , as long as Vq is sufficiently small (in comparison to an appro- 
priate energy scale contained in H, such as Pq in the case of the Klein-Gordon equation). 
We therefore expect it to have properties similar to Gp, as described in Section ID. 

However, Vq is not set to zero exactly and this in fact means that Gy has two properties 
not possessed by Gp- First, the non-zero Vq produces a suppression for widely separated 
initial and final points. In a path integral representation of Gv{q,q'), the sum will be 
dominated by classical paths from q to q', to which one may associate the total parameter 
time r of the path. From the integral representation Eq.(4.6) it can be seen that Gy will have 
an overall exponential suppression factor of the form exp(— Vor). Thus, propagation with 
Gv will suppress configurations q and q' connected by classical trajectories of parameter time 
duration of greater than I/Vq. Second, recalling that the class operator is closely related 
to the (S-matrix for scattering off a complex potential there will be some reflection involved 
and this appears in properites of Gy - it will not exactly possess the Feynman properties of 
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propagating positive frequencies to the future negative frequencies to the past, but may, for 
example, propagate some positive frequencies to the past. However, this "non-Feynman" 
propagation will be small for sufficiently small Vq, as we will see in the next section, so this 
possibility will be ignored. 

These properties of Gy are crucial to understanding the properties of Ca. First note that 

{q\PGvP\^x) = -^ ^ (t-^q! Mq)Gv{q, q')K'^x{q') (4.14) 

In this expression the propagator Gv{q, q') propagates from the boundary E to a point q on 
the interior. The fact that Gy is, approximately, a propagator of the Fcynman type means 
two things. First, there will be initial and final points q, q' for which the propagator is 
very small (due to "propagation outside the lightcone", discussed in Section ID). Secondly, 
when it is not small, one would expect it to involve only ingoing modes at the boundary (to 
the extent that reflection is ignored). Hence the PGy terms effectively restrict the current 
operator P at the boundary to ingoing modes only - exactly the result we are looking for. 
We therefore replace P, with the current operator (P)i„ involving ingoing modes only 

Ca|*a> = 2W^5{H - A) PGv{P)in\'^x) (4.15) 

Now note that since the PGy terms have done their job of selecting the ingoing modes, 
they may be eliminated. More precisely, we write P — 1 — P , where P is the projector onto 
the region outside A, and noting that that 

Vq5{H - \)Gv = 5{H - A) (4.16) 

we find 

Ca|*a) = 27r(^(i/-A)(P),„|*A) 

- 2t:5{H - X)PGviP)i„\^x) (4.17) 

Consider the second term on the right-hand side. It consists of ingoing modes at the bound- 
ary of A propagated with Gy to a final point which is outside A. Semiclassically, this 
corresponds to paths which have to traverse the entire width of A, so if A is sufficiently 
large, this will take a long parameter time r and, as discussed above, Gy will contain a 
suppression factor of exp(— Vqt) in comparison to the first term in Eq.(4.17). This is exactly 
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analogous to the suppression of the second term in the classical case, Eq.(4.3). We thus 
deduce that 

CaI^a) ~ 27r5(i/ - A)(P),„|^a) (4.18) 

Note that the approximation leading to dropping the second term in Eq.(4.17) also ensures 
that the result is independent of Vq, like the classical case in the appropriate hmit. Eq.(4.18) 
is the main result of this section. 

Eq.(4.18) may also be expressed in terms of {P)out defined in Eq.(4.12). To see this, note 
that, since P = i[H, P], we have 

5{H - A')P|*a) = ^(A' - X)S{H - A')P|^a) (4.19) 

This is zero for A = A', as long as (A|P|A) is well-defined, where |A) are eigenstates of H. (It 
may not be well-defined if P projects onto an infinite region - see Appendix A). Eq.(4.12) 
then implies 

S{H - A)(P),„|^a) = S{H - X){P)out\'^x) (4.20) 

so Eq.(4.18) may be expressed in terms of either the ingoing or outgoing flux or both. 
It is also useful to note that the right-hand side of Eq.(4.18) may be written 

27r5{H - A)(P),„|*a) = /e|*a) (4.21) 

where 

/oo 
dt e'"\P),ne-'"' (4.22) 
-oo 

This is an intersection number operator for ingoing flux at the boundary E of A. Hence 
the class operator is essentially the intersection number /e, and clearly commutes with H. 
It is classically equal to 1 for trajectories with ingoing flux at E (with no more than two 
intersections of the boundary, in the approximation we are using), and zero for trajectories 
not intersecting E. Of course, one may have guessed this approximate formula for the class 
operator, but a class operator is fundamentally defined as a product of projectors (or quasi- 
projectors) and the derivation given here makes it clear how this obvious guess arises from 
the fundamental definition. Note also the the class operator is hermitian in this case, so 
there is no need to consider the modified propagator Eq.(3.11). 

The form Eq.(4.18) may be used to check for decoherence of histories in specific models 
and we will see this later. When there is decoherence of histories, the probabilities are given 
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by the average of a single class operator, Eq.(1.19), which in this case reads 

(*a'|Ca|*a) = 27r(*A|(P)m|*A) ^(A - A') (4.23) 

Following the induced inner product prescription, we drop the 5-function on the right and 
then set A = A' = 0, so the probability in terms of a solution l^*) of the Wheeler-DeWitt 
equation is 

(*|CA|*)pfe,. = 27r(*|(P),„|*) (4.24) 

This is essentially the ingoing Klein-Gordon flux on the boundary of A, as expected. (The 
factor of 27r relates to the induced inner product as described in Appendix A). 

A more detailed calculation in a specific model is required to see that the above argument 
works. In particular, it is necessary to show that the various requirements on Vq can be 
simultaneously met - it has to be small enough for Gy to function as a Feynman propagator 
and for reflection to be neglected, but large enough to ensure independence of Vq and the 
dropping of the second term in Eq.(4.17). We will see in the models below that this is indeed 
possible. 

B. Two Regions 

We flrst consider the class operator Eq.(3.29) for two regions in the classical case. Suppose 
a classical trajectory enters Ai at ta, leaves at tb, enters A2 at tc and leaves at td (with no 
more than two crossings of the boundaries of each region). See Figure 5. Then we have 

p2{t')Pi{s') = [5{t' - g - S{t' - U)] [5{s' - ta) - 5{s' - h)] (4.25) 

It is easy to see that only one of the two terms in Eq.(3.29) contributes - the other term 
corresponds to the time-reversed trajectory. One can also see that the 5-functions at t' — td 
and s' — ta make no contribution to the integral. Hence the only contribution comes from the 
time of departure tb from Ai and the time tc of entering A2, clearly sufficient to characterize 
histories which spend time in both regions without overcounting. Evaluating the integral 
we find 

= (1 - e-^^"^*'*-*^)) (1 - e-^^"^*''-*'")) (4.26) 

which is approximately 1, as required, as long as the trajectory spends sufficiently long in 
each region. For trajectories entering only one region, the P term associated with the other 
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region is zero, so the class operator is zero. Classically, the class operator therefore has all 
the desired properties. 

In the quantum case, we follow steps similar to the case for one region considered above. 

We consider the action of the class operator CaiA2) defined in Eq.(3.29) on an eigenstate 
I^'a) of H. Changing variables from s, s\ t' to 

s — t — s 
s' ^ s'-s 

f ^ t-t' (4.27) 
the t integral may be done with the result 

/•OO PS ps — s' 

CaiA^I^a)^ - 27r / ds ds' / di' exp (-1/20? - lios') 

JQ Jo Jo 

+ (1^2) (4.28) 

Now note that 

^ -T4oe-'^^"'(Pi)o„t|*A) (4.29) 

The first approximation arises because the combination Pi exp{—iHs')Pi is propagation from 
inside Ai to the boundary with only positive s', which means that there are only outgoing 

modes at the boundary. The second approximation arises because, once the restriction is 
made to outgoing modes, the projection Pi may be dropped (similar to dropping the second 
term in Eq.(4.17)). Similar approximations may be made with the terms 5{H — X)V2e~'^^^' P2. 
We thus obtain 

POO PS ps — s' 

CaiA^I^a)^ 27rl/iol/2o / ds ds' exp (-I/20? - V^ios') 

Jo Jo Jo 

X S{H - X){P2)ine-'^''-^^^'-''-''\Pi)out\^x) 

+ (1^2) (4.30) 
The integrals may be carried out with the result, 

CaiA.I^a) = ^TrSiH - A) [{P2)inGF{X){Pl)out + (1 ^ 2)] I^T;,) (4.31) 
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where G'f(A) is the Feynman propagator at eigenvalue A, 

/•oo 

Gf{\) = / dt e-'^^-^-''^ 
Jo 

We may again attempt to interchange ingoing and outgoing flux operators, along the lines 
of Eq. (4.20). However, since 

{H-X)GF{X) = -i (4.33) 

there are extra terms of the form (P2)m-Pi and (Pi)m-f2, but these are zero since they are 
products of operators locahzed about spatially distinct regions. We thus find 

CaiA.I^a) = 27rS{H - A) [{P2)inGF{X){Pi)in + (1 ^ 2)] I^a) (4.34) 

(or equivalently in terms of outgoing flux). 

At this stage it is of interest to examine the difference between this propagator and the 
modified version, Eq.(3.11), invariant under reversals of parameter time, since the presence 
of Gf(A) in Eq.(4.34) rather than S{H — A) means that it is not invariant. It is easy to see 
that the modification Eq.(3.11) effectively produces the replacement 

Gf{X) - I {Gf{X) + GUA)) = l{27:)d{H - A) (4.35) 
and we obtain for the modified propagator 

C'a.aJ'^x) = ^{2nfS{H - X){P2)inS{H - X){P,U<ilx) 

+ ^{27rf 6{H - X){PiUS{H - XmU^x) (4.36) 



This result may also be expressed in terms of intersection number operators, 

1 
2 

It is now easy to guess the form of the class operator for n regions. It is. 



Ca,a. = o(We,+/e,/e.) (4.37) 



C'n^ —, • • • -^s„ + permutations) (4.38) 

77/. 

where "permutations" means add all possible permutations of the n regions, to give a total 
of n\ terms. This final result has a particularly natural form which also suggests that it 
may have simple path integral representations. This will be explored elsewhere. (See also 
Ref.[41] for similar expressions, derived using a detector model for quantum cosmology). 
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V. SOME SIMPLE EXAMPLES 



To verify the properties of the class operator Ca outhned above, we compute it exphcitly 
in some simple examples. We first consider a one-dimensional example. 

A. A One-Dimensional Example 

We take the Hamiltonian to be the free particle Hamiltonian H = /2m and the region 
A to be the interval [— L, L] . We take the initial state to be an energy eigenstate 

<-lf> = p^^"' (5-1) 
with p > and energy E — /2m. (In this case, unlike the Wheeler- De Witt and Klein- 
Gordon equations, the eigenvalue of H is the physical energy, so we do not take = at the 
end of the calculation. However, the inner products are still regularized by taking different 
values of E on either side of any inner product.) 

Our aim is to confirm in this example the properties of the class operator Eq.(4.5), which 
in this case we write as, 

CaIp) = 2W^5{H - E)PGvP\p) (5.2) 

P is the projection operator onto the region [— L, L] and may also be represented in terms 
of the region's characteristic function fiix), so P — fiix). It follows that 

p = ^ \p5(x + L) + 5ix + L)P - pS(x - L) - 5(x - L)p] (5.3) 
2m 

For the p > initial state considered here, we identify the first two terms, at 

{P)in and the last two terms, at a: = L, as —{P)out^ 

P = (P),„ - {P)out (5.4) 

We may now write 

{x\PGvP\p) = J dyGv{x,y){y\P\p) 

=/i(a;) f[Gy(x,y)^eM -[Gy(x,y)^eH ] (5.5) 



2mV27r 

The Green function may be calculated explicitly using the form Eq.(4.6), which may be 
written 

Gv{x,y) = J^ dt e^plim^^f^ + i{E + iVo)t] (5.6) 
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The integral may be carried out explicitly [42] , with the result 

Gv{x, y)^^ [9{x - t/)e^(^-^)« + 9{y - x)e-^(^-^)«] (5.7) 
where we have introduced the complex momentum 

Q = ^/2m{E + iVo) (5.8) 

The interesting case is that in which Vq <^ E (which is the condition for negligible reflection 

in the scattering off the complex potential) and we can then expand 

1/2 



Q = V2mE (^1 + ^^^ 



= p + + (5.9) 
P 

where p — y/2mE. This means that the Green function has the form 

Gv(x,y) ^ - \e(x - y)e'^''-y^P + 0(y - x)e-^(^-^)fl exp ( -^\x - y\) (5.10) 
P \ P J 

It therefore has the form of Feynman type Green function, but there is also exponential 

suppression for large separations of initial and final points, as anticipated. 

To evaluate Eq.(5.5), we consider the action of the Green function Gy on a plane wave 

state \p) with p > 0. We have 

^Gy{x,y)"dl(^^ = e{x-y)^^^e^v{ipx + i{x-y){Q-p)) 

(Q _p) 

- 0{y - x) exp {-ipx + 2ipy + i{y - x){Q - p)) (5.11) 

For E Vq, this has the approximate form, 

—Gv(x,y)d^e'Py = Oix - y) e^^^' exp ( -"^W ~ y\\ 
2m \ P J 

- Ai^e{y - x) e-^f^+2^f2/ exp (^-"^k - (5-12) 

The first term corresponds to the familiar property of consisting of positive momentum 
ior X > y (with the exponential suppression factor mentioned above). The second term 

is not a familiar property of the Feynman Green function since it corresponds to negative 
momentum for x < y (which we do not expect since the incoming state consists only of 
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positive momentum). This actually corresponds to reflection, and this term will be small 
for Vo <^ E, which we assume. 

These results show that the second term in Eq.(5.5), the term corresponding to {P)out 
may be dropped, as expected. We have now confirmed that 

CaIp) ~ 2t:V^5{H - E)PGv{PUp) (5-13) 
Following the general discussion of Section 4, we write P — 1 — P and obtain 

Ca\p) ^ 2t:5{H - E){P)M 

- 27rVoS{H - E)PGv{P)in\p) (5.14) 

We need to show that the second term is small. We have 

{x\PGv{P)in\p) = —^hix) [Gy(x,y)^eH (5.15) 



y=-L 



2mV27r 

where fiix) is the characteristic function for the region outside A, so x > L and x < —L. 
The properties of Gy deduce above confirm that this expression is small. For x > L, 
there is exponential suppression of the Green function since |x — y| > 2L. The exponential 
suppression is of order exp(— Vqt), where 



r = ^ (5.16) 
P 

which is the time for a classical trajectory to traverse A. On the other hand, values of 
X < —L can only be reached by reflection, and this is small since Vq <^ E. 

We have therefore confirmed the desired result that the class operator has the form 

GA\p)^27r5{H-E){PUp) (5.17) 

and this result is valid in the regime 

-<^Vo^E (5.18) 

T 

(This regime is a commonly encountered one in studies of the arrival time in non-relativistic 
quantum mechanics [29, 31]). The outer parts of the inequality imply that Er ^ 1, or 
equivalently, 

pL > 1 (5.19) 
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which is easily satisfied for sufficiently large L and p. It is essentially the condition that the 
size of A is much greater than the wavelength of the quantum state. 
Since 

2tx5{H - E) =Gf + G\ (5.20) 
(where is the Feynman Green function, obtained from Gy by letter Vq — > 0+), we have 

ix\5iH - E)\y) = ^ Uvi^-y) + ^-ip{^-y)\ (5.21) 
p 

and it follows that Eq.(5.17) may be evaluated with the trivial result 

CaIp) « \p) (5.22) 

This is completely expected since every plane wave enters the region A. 

However, for the class operator operating on the initial state Eq.(5.17) to have a more 
interesting and non-trivial form, we need to look at a two-dimensional example. This will 
be useful for the purposes of examining the decoherence properties of some simple states. 



B. A Two-Dimensional Example 

We again take a free particle but in two dimensions, with Hamiltonian H — (p^-|-p|)/2m, 
and we take the region A to be a rectangle of sides 2Li and 2L2 centered at the origin. We 
use coordinates xi,X2 and the characteristic function of A is 

fA{,Xi,X2) = fLA^l)fL^{x2) (5.23) 

We take the initial state 1^^) to be a plane wave in the 1-direction, e'^^^^ with pi > 0, so 
there is no momentum in the 2-direction. The analysis of the one-dimensional example may 
be used to show, very easily, that the class operator is again of the form Eq.(5.17), with 

{Phn = ^ \pi5{x, + Li) + 6{x^ + L,)p,] fL,{x2) (5.24) 

Denoting the matrix elements of S{H — E) by G'(xly), we have 



(x|Ca|*) = 2T:{xr,X2\5{H - E){P),^\^) 
~ 2m 



dy2 \G{x^,X2\yi,y2)K^'''''^'] (5-25) 

J-L2 '- J2/i=-Li 
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Note that, due to the suppression effect in Gy, the spatial size 2Li of A in the Xi-direction 
has essentially dropped out of the expression for the class operator and it is expressed now 
only in terms of the incoming flux at Xi = —L. 

Eq.(5.25) may be approximately calculated, at some length, but the form of the result for 
Ca|^) is easy to anticipate on general physical grounds. We still have the result Eq.(5.20), so 
we first consider the action of Gp attached to the initial state as in Eq.(5.25). The physical 
situation therefore consists of an incoming plane wave exp(ipiXi) from xi < encountering 
a gate of width 2L2 in the X2 direction located at Xi — —L. This situation is very similar to 
the Mott analysis of the particle detection in a cloud chamber [41, 43]. The localization by 
the gate to spatial width 2L2 produces a momentum uncertainty in the 2-direction Ap2 of 
order I/L2. Since the momentum in the 1-dircction is pi, it is easy to see that this produces 
a spreading in xi > of angular size of order 1 / (P1L2) . This will be very small for sufficiently 
large L2. If we now consider also the action of G^p on the same initial state, it produces 
essentially the same effect in xi < 0. 

The final result is therefore that Ca|^) consists of a plane wave exp(ipiXi) locahzed on a 
tube of width 2L2 around the Xi-axis, with the tube opening out by a small angle l/{piLi) 
as Xi — > ±00. See Figure 6. To the extent that the angular spreading can be ignored, the 
result has the very approximate form 

{xi,X2\CA\^)^e'P^^^fL^{x2) (5.26) 

where /^^ ^ function localized to a width L2 around the origin (not necessarily the same 
as the exact window function used above). 

This approximate result may be used to address the decoherence properties of the initial 
state 1^) defined above. The class operator for not entering A is clearly 

(xi,X2|CA|^)~e*^i^7L(^2) (5.27) 

where /^^ = 1 ~ /l2' C'a|^') is localized outside the tube defined in Eq.(5.26). We 
immediately see that 

(^ICaCaI*) ~ (5.28) 

quite simply because the two class operators are spatially localized about different non- 
overlapping regions. There is therefore approximate decoherence of histories in this case. 
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Another interesting case is to take the same region A but take a superposition of initial 
states 

|vl>) = i=(|vI/,) + |xl;2)) (5.29) 

where |^i) is a plane wave in the 1-direction, exp(ipia;i), considered above, and 1^2) is a 
plane wave in the 2-direction, exp{ip2X2)- Each state individually produces decoherence of 
histories, as shown above, but there will be cross terms in the decoherence functional, 

(^ICaCaI*) ~ ^(*2|C'aCa|*i) + ^(*i|C'aCa|*2) (5.30) 

The first term is approximately, 

(^'2|C'aCa|^'i) ~ J dxidx2 fL^{x2)fLiixi)exp{ipiXi+ip2X2) (5.31) 

It is easy to see that this expression will be very small as long as piLi » 1 and P2L2 ^ 1, 
due to averaging out of oscillations. There will therefore be decoherence of superposition 
states as long as the region A is sufficiently large. 

The essential reason for decoherence in this case is that the orthogonality of |^i) and 
1^2) is little disturbed by the action of the class operators as long as A is large enough. 
(Recall that interference effects in the double slit experiment will only arise for narrow slits 
very close together). For regions A of small size, it may be necessary to consider coupling 
to an environment to produce decoherence of histories. This will not be explored here. 



C. The Klein-Gordon Equation and the Probability for Surface Crossing 

We now consider the problem of assigning probabilities to histories of the relativistic 
particle in flat space which cross or do not cross the spacelike surface = in the spatial 
region E. This is shghtly different to the cases considered so far since we need to show exactly 
how the results derived above for entering regions A may be used to derive probabilities for 
crossing surfaces. This example also involves the effects of the signature of the metric. A 
decoherent histories analysis of the Klein-Gordon system was given in detail in Ref.[10], 
but it encountered difficulties in definition of the class operators. A somewhat heuristic 
definition was used, with satisfactory physical results, but the underlying origin of this 
definition was not clear. Here, we show how this heuristic definition arises from the more 
rigorous definitions used in this paper. 
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To define a class operator for crossing the spacelike surface E , we regard the surface as 
part of the boundary of a spacetime region A. For simplicity we take the three-dimensional 
region E to be the interior of a 2-sphere at — 0. We then define A by requiring that it is a 
null cone in > with S its base. (A two-dimensional version of this is depicted in Figure 
7). The class operator is again given by Eq.(4.18), with |^) a solution to the Klein-Gordon 
equation which may contain positive and negative frequencies, 

^ = ^+ + ^- (5.32) 
The flux operator {P)in consists of contributions from the two parts of the boundary of A, 

inn = {P)l + {P)T' (5.33) 
The flux operator on E is clearly just the positive frequency flux at this boundary 

{P)l = {e{po)m^') + 5{xy{po)po) Mx') (5.34) 

where fT.{x^) is a window function localized in E. The remaining part of the flux operator 
describes the incoming flux on the null cone. Classically, since this surfaces is null, every 
incoming trajectory crossing the cone will cross E. An analogous result is approximately 
true in the quantum case, since propagation outside the lightcone is suppressed. Hence, 
these flux operators are approximately equal to the negative frequency flux at E, 



{P)tr ~ {0{-Po)PoS{x') + S{xy{-po)po) h{x') (5.35) 
The class operator involves contributions from boths of these terms and may be written 

Cj:\^x) = 27iS{H - A) (|^)o|<5(£°) + S{x')\po\) /e(£OI*a) (5.36) 
Noting that {x\5{H)\y) is the same as the the Klein-Gordon propagator G*^^) = -|- G~ 



(where are the positive and negative frequency Wightman functions [10]), we have, 
setting A = 0, 

(x|Cs|*) = 2m J £y 

This agrees with the deflnition in Ref.[10]. It consists of the initial state attached to G*^^) 
with a "sign- adjusted" Klein-Gordon inner product to take care of the negative frequencies. 
It was shown in Ref.[10] that there is approximate decoherence between histories which cross 
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GW(x,y)ao*+(y)-G«(x,y)^*-(y) 



(5.37) 



x° = in E or outside E, again because propagation outside the lightcone is suppressed. 
The resulting probabihty for crossing is 



This is the expected result and is positive. 

Hence, the decoherent histories approach described in this paper yields the expected 
results for the Klein-Gordon equation, improving on the heuristic analysis of Ref.[10]. Note 
that the method used above to relate probabilities for surface crossings to probabilities for 
entering regions is specific to the dynamics of the rclativistic particle in fiat space. More 
complicated systems will require a modified approach adapted to their particular dynamics. 

VI. WKB REGIME 

The most important case in which to check the ideas developed above is in the WKB 
regime. In the oscillatory regime, the solutions to the Wheeler-DeWitt equation have the 
form 



where R and S obey Eqs.(1.24), (1.25) as described earlier. More generally, the wave function 
is a superposition of WKB wave functions but we consider first the case of a single term. The 
heuristic interpretation of such states was described in Section ID. Our aim is to show that 
the decoherent histories analysis reproduces the heuristic scheme. WKB states are locally 
plane wave states of the type considered in the previous Section, so we will appeal to that 
analysis to understand the properties WKB states. 

We consider the action of the class operator for a single region A on a WKB state. We 
first consider the class operator Eq.(4.18) acting on WKB state (regularized by making it 
an eigenstate of H with eigenvalue A) 




(5.38) 



* = Re 



.iS 



(6.1) 



(?|Ca|*a) = 2n{q\S{H-X){PU^^) 




(6.2) 



Here Ej„ denotes the sections of the boundary where the fiux is ingoing. 



n-VS <0 



(6.3) 
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where n" is the outward pointing normal. The key property of the WKB wave function 
^\{q') is that it has momentum p — 'S7S{q') at each point q' on the boundary. When the 
operator 5{H — A) is apphed, from the representation Eq.(4.8), we see that it's effect is 
to evolve the state "^xiq') (restricted to forwards and backwards in parameter time, 
and then integrates over all times. Semiclassically, the evolution of the state \E'a(q'') will 
be concentrated along the classical trajectories defined by initial positions q' in Ej„ and 
momenta p — VS{q'), with some spreading of the wave packet, but this will be small if A is 
reasonably large. (This is analogous to the model of Section 5B). 

We thus see the following: the wave function ((/ICaI^a, ) in Eq.(6.2) is spatially localized 
around the tube of classical trajectories passing through A with momenta p = VS (depicted 
in Figure 2). This may be approximately written in the alternative form 

(g1CA|*A)~^^(T-A-e)i?e*^ (6.4) 

where e > is a small parameter to regularize the ^-function at zero argument. Here, 
T/\{q) is the parameter time spent by the classical trajectory qci{t) (with initial value q and 
momentum p = VS{q)) in the region A and may be written 

/oo 
dt fAiqciit)) (6.5) 
-oo 

This has the property that 

• Vta = (6.6) 

since VS ■ V simply translates along the classical trajectories. It follows from Eq.(1.25) that 
Eq.(6.4) is in fact a WKB solution to the Wheeler-DeWitt equation, since the ^-function 
essentially modifies the prefactor R but in a way that it still satisfies Eq.(1.25). (Of course, 
this is expected because Ca commutes with H). 

Eq.(6.4) is a very useful result and allows us to check for decoherence very easily. The 
action of the class operator for not entering A is clearly 

(^IC-aI^a) -^(e-TA) i?e^^ (6.7) 

which is a WKB state locahzed on the set of trajectories not entering A. It immediately 
follows that 

(*a'|C'aCa|*a) ~0 (6.8) 
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since the two states Eqs.(6.4), (6.7) are localized about complementary regions. There is 
therefore approximate decoherence of histories for a single WKB packet and for histories 
entering or not entering a single region A, as long as A is sufficiently large. 

The key reason for the decoherence with a single WKB packet is related to the approx- 
imate determinism of the WKB wave functions: fixing values of position to lie on Ej„ also 
fixes the momenta, since p — WS, so that Eq.(6.2) is concentrated along a tube of classical 
trajectories. 

More generally, the initial state will be a superposition of WKB wave packets. 



The component states in this sum are typically approximately orthogonal to each other as 
long as the phases Sk are sufficiently different. (This will depend on the detailed dynamics of 
the model). Following the analogous example in the simple models of Section 5B, we would 
expect that the class operators will not disturb the approximate orthogonality of these 
states as long as the region A is sufficiently large - again we expect the cross terms in the 
decoherence functional to average to zero because of oscillations, as in Eq.(5.30). Therefore, 
superpositions may in practice be treated as mixtures at sufficiently coarse-grained scales. 
Note that this statement also applies the special state. 



which arises from the Hartle-Hawking "no boundary" proposal [44]. That is, the interference 
between the two terms may be neglected. 

Given decoherence, the probabihties are given by the general expression Eq.(4.24). It 
then follows that the probabilities for entering A coincide with Eq.(1.29), the sought-after 
result. 

Now consider the case of probabilities for histories entering two regions, as described by 
the (modified) class operator Eq.(4.36). Since the two- region class operator is a sum of 
products of one-region class operators, its effect on the WKB wave functions is easy to see. 
The action of a single class operator gives Eq.(6.4). But since this is still a wave function of 
the WKB type, the action of a second class operator yields 




(6.9) 



k 



(6.10) 




iS 



(6.11) 
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That is, it is a WKB wave function but restricted in such as way that it's flux passes through 
both regions. See Figure 8. (It is easy to see that in this case the action of the unmodifled 
class operator Eq.(4.34) is the same.) It is again easy to see that there is decoherence of 
histories and the probabihty is given by an expression of the form Eq.(1.29), but with the 
integral over the subset of incoming flux at Ai which goes on to intersect A2. 

Prom the above observations one can also see why the exponential WKB wave functions 
will not lead to decoherence of histories. The exponential wave functions have the form 

* = e"^ (6.12) 

where / is real. The key difference between states of this type and the oscillatory type 
Eq.(6.1) is that they do not have a correlation between positions and momenta [28]. One 
would therefore expect the evolution of the state Eq.(6.2) to be spread all over the conflgu- 
ration space and not concentrated around a particular region. That is, these states do not 
have the approximate determinism of the oscillatory states. The states Ca|^) and Ca\^) 
would then not be approximately orthogonal so there will be no decoherence of histories. 

The decoherence of histories described here has arisen because of the approximate deter- 
minism of the oscillatory WKB states, together with the approximate orthogonality proper- 
ties that arise when the regions A are sufficiently large. At finer grained scales, decoherence 
of histories may only be possible in more complicated models in which there is an envi- 
ronment of some sort. Models along these lines, in more basic approaches to quantum 
cosmology, have been considered previously [45], and one might expect that they may be 
adapted to the decoherent histories approach to quantum cosmology. (See also Ref.[ll]). 

In summary, we have derived from the decoherent histories approach the probabilities 
normally used in the heuristic WKB interpretation. To address these issues in more detail 
will require more specific quantum cosmological models. This will be considered elsewhere. 

VII. DISCUSSION AND FURTHER ISSUES 

We have presented a properly defined quantization procedure for quantum cosmology us- 
ing the decoherent histories approach to quantum theory and derived from this the frequently 
used but heuristic WKB interpretation, involving fiuxes of the WKB wave function. 

The key idea was to use a complex potential to define the class operators for not entering a 
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region of configuration space. Tfiis method is adequately justified by its successful use in the 
arrival time problem in non-relativistic quantum theory. We showed that the class operators 
defined in this way have all the desired properties - they have the correct classical limit, are 
compatible with the constraint equation, and do not have difficulties with the Zeno effect. In 
a semiclassical approximation, they have an appealing form in terms of intersection number 
operators. They give sensible results in simple models and there is approximate decoherence 
of histories for certain types of initial state at sufficiently coarse grained scales. 

Future papers will address the more detailed application of this approach to specific 
models and will also undertake a comparison of the decoherent histories approach described 
here to other approaches [5]. 
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APPENDIX A: OPERATORS AND INNER PRODUCTS FOR THE KLEIN- 
GORDON EQUATION 

In this appendix we outline the relationship between the induced inner product and Klein- 
Gordon inner product. This has been described elsewhere (see, for example, Ref.[9]), but 
the derivation given here is different and will help give some insight into the form of some 
of the results of Section 4. 

We use = (a;°, x) to denote spatial coordinates and take the signature of spacetime to 
be ( — h ++). The Hamiltonian for the Klein-Gordon system is 



The induced inner product begins by looking at eigenfunctions of H with eigenvalue A, which 
we write as 



(Al) 



(A2) 
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where po — a/p^ + m'^ — A and s = ±1 labels the solutions as positive or negative frequency. 
We work with three different inner products. The first is the standard Schrodinger inner 
product, 

(Ap^lA'pV) = J d'x ^l^,{x)^y^,s'{x) (A3) 
and the eigenfunctions are normalized with this inner product according to 

(Aps|A'p's') = 5{X - A') 5{p - p') 5ss' (A4) 

The goal is to relate this to the Klein- Gordon inner product Eq.(1.6). One way to do this 
is to integrate out x° in Eq.(A3), but here we proceed differently. 

We introduce the projection operator P — 9{x^) onto x° > and consider the current 
operator 

/oo 
dt P{t) (A5) 
-oo 

On the one hand, we can carry out the integral with the result 



e{f) (A6) 



where e(p°) is the signum function operator. On the other hand, we can differentiate P with 
the result 

P{t) = i[H,9{x^)] 

= + (A7) 

so that 

/oo 
dt e^^* + e-^^* (A8) 

-oo 

We now insert the identity 

e(p°)J=l (A9) 
with J given by Eq.(A8) in the inner product Eq.(A4), with the result 

(Aps|A'pV) = (Aps| e{f)J |A'pV) 

= 27r5(A - A') s (Aps| {p%{SP) + |Ap's') (AlO) 
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We now identify the inner product expression on the right-hand side as the Klein-Gordon 
inner product between states with the same A, 

(*Ap.,*Apv)i^G = (Aps|(/(5(£°) + 5(£V)|ApV) 

^ t [ d'x **p,(x)^*Ap'.'(2:) (All) 

(where = —pQ = ido). Therefore the relationship between the Schrodinger inner product 
and Klein-Gordon inner product is 

(Aps| A'p's') = 27r5(A - A') s (*Aps, '^Xp's')kg (A12) 

Finally, the induced (or physical) inner product on eigenstates with the same A, is defined 
as outlined in Section 1, essentially by factoring out the (5-function 

(Aps|A'p's') = 6{X - A')(Aps|Ap's%,,j,, (A13) 

so we deduce that 

(Aps|ApV)pftj,, = 27rs (*Ap., *Ap'.')kg (A14) 

and we may set A = to obtain an inner product on solutions to the Klein- Gordon equation. 
The Klein-Gordon inner product is of course positive on positive frequency (s = 1) solutions 
and negative on negative frequency (s = —1) solutions, but the form of the above result, 
with the overall factor of s, ensures that the induced inner product is always positive. 

Note that some caution is required in relation to Eqs.(A7) and (AlO), since Eq.(A7) 
implies that 

(Aps|(/5(£°) + 5(x°)/) lA'pV) = i(Aps|[ii',^(x°)]|A'pV) 

= i(A-A')(Aps|e(x°)|A'pV) (A15) 

This appears to be zero when A = A' suggesting that Eq.(AlO) is zero. This is not in fact 
the case since the inner product expression on the right-hand side of Eq.(Al5) is actually 
divergent as A — > A'. 



[1] J.J.Halliwell, in Quantum Cosmology and Baby Universes, edited by S.Coleman, J.B.Hartle, 
T.Piran and S.Weinberg (World Scientific, Singapore, 1991) also available as the e-print gr- 
qc/0909.2566. These introductory lectures on quantum cosmology contain a detailed account 
of the WKB interpretation. 

42 



[2] C.Kiefer and B.Sandhoefer, gr-qc/0804.0672; D.L.Wiltshire, gr-qc/0101003. 
[3] J.B.Hartle, S.W.Hawking and T.Hertog, Phys.Rev.D77, 123537 (2008); Phys.Rev.Lett.lOO, 
201301 (2008). 

[4] B.DeWitt, in Gravitation: An Introduction to Current Research, edited by L.Witten (John 

Wiley and Sons, New York, 1962). 
[5] CRoveUi, Class.Quant.Grav. 8, 297 (1991); 8, 317 (1991); Phys.Rev. 42, 2638 (1990); 

Phys.Rev. 43, 442 (1991); D.Marolf, Class.Quant.Grav. 12, 1199 (1995); 12, 1441 (1995); 

P.Hajicek, J.Math.Phys. 36, 4612 (1995); M.Montesinos,C.Roveni and T.Thiemann, 

Phys.Rev. D60, 044009 (1999); M.Montcsinos, Gen.Relativ.Gravit. 33, 1 (2001); R.Gambini 

and R.A.Porto, Phys.Rev. D63, 105014 (2001). 
[6] K.Kuchar, gr-qc/9304012 and in Proceedings of the 4th Canadian Conference on General 

Relativity and Relativistic Astrophysics edited by G.Kunstatter, D.Vincent and J.Williams 

(World Scientific, Singapore, 1992); C.Isham, gr-qc/9210011. 
[7] S.W.Hawking, Nucl.Phys. B239, 257 (1984). 

[8] J.B.Hartle, in Proceedings of the 1992 Les Houches School, Gravity and its Quantizations, 
edited by B.Julia and J.Zinn-Justin (Elsevier Science B.V., 1995) 

[9] J.B.Hartle and D.Marolf, Phys.Rev. D56, 6247 (1997). 
[10] J.J.HalliwcU and J.Thorwart, Phys.Rev. D64, 124018 (2001). 
[11] J.J.HalliweU and J.Thorwart, Phys.Rev. D65, 104009 (2002) 
[12] J.J.HalliweU and P.Wallden, Phys. Rev. D73, 024011 (2006) 
[13] D.Craig and J.B.Hartle, Phys. Rev. D69, 123525 (2004). 
[14] J.Whelan, Phys.Rev. D50, 6344 (1994). 

[15] C.Anastopoulos and K.Savvidou, Class.Quant.Grav. 17, 2463 (2000); 22, 1841 (2005). 

[16] A.Ashtekar, J.Lewandowski, D.Marolf, J.Mourao and T.Thiemann, J.Math.Phys. 36, 
6456 (1995); A.Higuchi, Class.Quant.Grav. 8, 1983 (1991). D.Giulini and D.Marolf, 
Class.Quant.Grav. 16, 2489 (1999); Class.Quant.Grav. 16, 2479 (1999). F.Embacher, Hadronic 
J. 21, 337 (1998); N.Landsmann, J.Geom.Phys. 15, 285 (1995); D.Marolf, gr-qc/0011112. 

[17] M.Gell-Mann and J.B.Hartle, in Complexity, Entropy and the Physics of Information, SFI 
Studies in the Sciences of Complexity, Vol. VIII, W. Zurek (ed.) (Addison Wesley, Reading, 
1990) ; and in Proceedings of the Third International Symposium on the Foundations of Quan- 
tum Mechanics in the Light of New Technology, S. Kobayashi, H. Ezawa, Y. Murayama and 



43 



S. Nomura (eds.) (Physical Society of Japan, Tokyo, 1990). 
[18] M.Gell-Mann and J.B.Hartle, Phys.Rev. D47, 3345 (1993). 

[19] R.B.Griffiths, J.Stat.Phys. 36, 219 (1984); Phys.Rev.Lett. 70, 2201 (1993); Phys.Rev. A54, 

2759 (1996); A57, 1604 (1998). 
[20] R. Omnes, J.Stat.Phys. 53, 893 (1988). 53, 933 (1988); 53, 957 (1988); 57, 357 (1989); 

Ann.Phys. 201, 354 (1990); Rev.Mod.Phys. 64, 339 (1992). 
[21] R.Omnes, J. Math. Phys. 38, 697 (1997) 

[22] J.J.HaUiwell, in Fkmdamental Problems in Quantum Theory, edited by D.Greenberger and 
A.ZeiUnger, Annals of the New York Academy of Sciences, 775, 726 (1994). 

[23] J.J.HaUiwell, e-print quant-ph/0301117, to appear in Proceedings of the Conference, Deco- 
herence. Information, Complexity, Entropy (DICE), Piombino, Italy, September, 2002 (edited 
by T.Elze). 

[24] J.J.HaUiwell, Phys.Rev. D60, 105031 (1999). 

[25] C. Isham, J.Math.Phys. 23, 2157 (1994); C. Isham and N. Linden, J.Math.Phys. 35, 5452 

(1994); 36, 5392 (1995). 
[26] J.J.HaUiwell Nucl.Phys. B266, 228 (1986); B286, 729 (1987). 

[27] J.J.HaUiwell, Phys.Rev. D38, 2468 (1988); J.J.HaUiwell and J.B.Hartle, Phys.Rev. D41, 1815 

(1990); D43, 1170 (1991); C.Tcitclboim, Phys.Rev. D25, 3159 (1982); D28, 297 (1983). 
[28] J.J.HaUiwell, Phys. Rev. D36, 3626 (1987). 

[29] J.G.Muga, R.Sala Mayato and I.L.Egusquiza (eds). Time in Quantum Mechanics (Springer, 
BerUn, 2002); J.G.Muga and C.R.Leavens, Phys.Rep.338, 353 (2000). 

[30] J.G.Muga, S.Brouard, D.Macias, Ann. Phys. (NY) 240, 351 (1995); Ph.Blanchard and 
A.Jadczyk, Helv. Phys. Acta 69, 613 (1996); J.P.Palao, J.G.Muga, S.Brouard, A.Jadczyk, 
Phys. Lett. A233, 227 (1997); A.Ruschhaupt, J.Phys. A35, 10429 (2002); J.G. Muga, J.P. 
Palao, B. Navarro, I.L. Egusquiza, Phys. Rep. 395, 357 (2004); J.J.HaUiwell, Prog.Theor.Phys. 
102, 707 (1999). 

[31] J.J.HaUiwell and J.M.Yearsley, Phys. Rev. A79, 062101 (2009). 
[32] J.J.HaUiwell and Y.Yearsley, quant-ph/0903.1958. 

[33] N.Yamada and S.Takagi, Prog.Theor.Phys. 85, 985 (1991); 86, 599 (1991); 87, 77 (1992); N. 

Yamada, Sci. Rep. Tohoku Uni., Series 8, 12, 177 (1992); Phys.Rev. A54, 182 (1996). 
[34] J.J.HaUiwell and E.Zafiris, Phys.Rev. D57, 3351-3364 (1998). 

44 



[35] J.B.Hartle, Phys.Rev. D44, 3173 (1991). 

[36] J.Echanobe, A. del Campo and J.G.Muga, Phys. Rev. A77, 032112 (2008). 
[37] J.Yearsley, in preparation. 

[38] A.Aucrbach and S.Kivelson, Nucl.Phys. B257, 799 (1985); P. van Baal, in Lectures on Path 
Integration, edited by H.A. Cerdeira et al. (World Scientific, Singapore, 1993) [also available 
from the website littp://www. lorentz.leidenuniv.nl/vanbaal/HOME/publ.html]; 

[39] J.J.Halliwell and M.E.Ortiz, Phys.Rev. D48, 748 (1993). 

[40] J.J.Halliwell, Phys.Lett. A207, 23 (1995). 

[41] J.J.Halliwell, Phys.Rev. D64, 044008 (2001). 

[42] L.S.Schulman, Techniques and Applications of Path Integration (Wiley and Sons, New York, 
1981), chapter 7. 

[43] N.F.Mott, Proc.Roy.Soc A124, 375 (1929), reprinted in Quantum Theory and Measurement, 
edited by J.Wheeler and W.Zurek (Princeton University Press, Princeton, New Jersey, 1983). 
For further discussions of the Mott calculation see also, J.S.Bell, Speakable and Unspeakable in 
Quantum Mechanics (Cambridge University Press, Cambridge, 1987); A.A.Broyles, Phys.Rev. 
A48, 1055 (1993); M.Castagnino and R.Laura, Int.J.Theor.Phys. 39, 1737 (2000); J.Barbour, 
Class.Quant.Grav. 11, 2853 (1994); 11, 2875 (1994). 

[44] J.B.Hartle and S.W.Hawking, Phys. Rev. D28, 2960 (1983). 

[45] J.J.Halliwell, Phys.Rev. D39, 2912 (1989); C.Kiefer, Class.Quant.Grav. 4, 1369 (1987); 
T.Padmanabhan, Phys.Rev. D39, 2924 (1989). 



45 



Figure 1. Given a solution ^' to tlie Wheelcr-DeWitt equation, what is the probability that 
the system enters a series of regions Ai, A2 • • • in configuration space, or crosses a surface 
E, at any stage in the system's entire history? 




Figure 2. A WKB wave function with phase S corresponds to a set of classical trajectories 
with tangent vector V-S". The probabihty for entering a region A is the amount flux of the 
wave function intersecting A. It may be expressed either in terms of the ingoing flux across 
Sj„ or equally, in terms of the outgoing flux across T,out- 
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Figure 3. The quantum arrival time problem in non-relativistic quantum mechanics. Given 
an initial state localized entirely in a; > and consisting entirely of negative momenta, what 
is the probability that the particle crosses the origin during the time interval [^1,^2]? 
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Figure 4. A classical trajectory q"'{t) intersecting A enters at parameter time ta and leaves 
at parameter time t^. 
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Figure 5. A classical trajectory intersecting Ai and A2 enters Ai at ta, leaves at tb, 

enters A2 at tc and leaves at td- 
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Figure 6. The class operator Ca acting on the initial state |^) produces a state consisting 
of a plane wave exp(ipia;i) localized along a tube of width 2L2 centred around the Xi axis, 
with spreading through a small angle as Xi — > ±00. 
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Figure 7. The class operator for a relativistic particle to cross a surface E may be calculated 
by regarding E as part of the boundary of the spacetime region A, a section of null cone 
whose base is E. In two dimensions, depicted here, it is a triangle with base E whose other 
two sides are the null surfaces Ei, E2. 
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Figure 8. The class operator CaiA2 for two regions operating on a WKB wave function 
produces another WKB state locahzed around the flux passing through both regions. This 
is the same as the flux entering Ai across the surface E12. 
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